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It is standard procedure to fit an applicable isotherm equation to water vapor adsorption data 
using the method of least squares in arriving at a value for the surface area accessible to the water 
molecule. The least squares technique has been extended in the present investigation to determine, 
in addition and simultaneously, a "best value" for the zero-humidity sample weight of the material. 
The application is equally valid for desorption insofar as the zero-humidity weight is concerned, although 
the derived value for "surface area" from desorption data will be over-estimated in the general case 
because of hysteresis. There is no limitation on the range of humidities since the method is not re- 
stricted to the BET equation (i.e., between 0.1 and 0.3 r.h.). In fact, good agreement with the zero- 
humidity points measured experimentally has been obtained even from drying curves in which the 
relative humidity has been confined to the region above SO percent. An iterative method is employed 
in the calculations for which computer assistance is especially adaptable. Fortran IV programs are 
included in the appendix whose use requires no extensive computer experience. A fraction of a second 
in computer processing time is all that is required for each determination. 
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1. Introduction 

What constitutes the dry weight of a solid has never 
been completely resolved. Some materials like hydrox- 
ide gels of aluminum, silicon, iron, etc., continue to 
lose water when evacuated at ambient temperatures 
and more so on heating, yet a trace of water always 
seems to remain. Other materials readily give up their 
water, but suffer irreversible changes in the process. 
This leads to the adoption of a concept of water of 
constitution as distinguished from excess water. 

A novel approach to this problem of dry weight 
prediction is portrayed in this paper based on the 
assumption that adsorbed water is certainly not water 
of constitution, and, consequently, is excess. To 
establish the precise point of dryness, therefore, it is 
only necessary to determine the weight of the solid 
under conditions where the amount of physically 
adsorbed water would be zero. This is determined by 
fitting vapor adsorption theory to the experimental 
data by an unorthodox application of the least-squares 
technique. 

Since the development in 1938 of the Brunauer, 
Emmett, and Teller (BET) multimolecular adsorption 
free-surface equation [l] 1 , many other adsorption 
isotherm theories [2, 3, 4, 5 6, 7] have appeared in the 
literature. Most of the resultant isotherm equations, 



I indicate the literature references at the end of this paper. 



despite their differences, shared certain important 
attributes: 

(1) They were expressable in some linear form. 

(2) Each equation possessed physically significant 
parameters which could be evaluated either from 
the slope or intercept of the experimental plot 
or a combination of both. 

(3) These parameters bore a simple functional 
relationship with the specific surface of the 
adsorbent. 

For these reasons, any such theory would appear to 
qualify as a candidate for prediction of dry weight in 
the method to be described. 

At the onset of this study, no a priori preference 
was made of any one adsorption theory over the others, 
although, in practice, the BET equation is probably 
more widely used than all of the others combined. 
This reason alone would justify its inclusion among the 
theories selected. In addition to the BET, the other 
adsorption theories tested in this work were: the 
Polarization theory [6], Harkins and Jura (H-J) theory 
[3], and Polanyi's Potential theory [6]. 

It will be shown that the same experimental points 
that are normally employed using least-square methods 
to determine only the surface area of an adsorbent 
can, at the same time, be called upon to furnish an 
additional piece of information; namely, its dry weight. 
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2. Theoretical Development 

(1) The BET [1] free-surface equation in its most 
familiar form may be represented as a linear function 
of the relative pressure: 



y = a + bx 



(i) 



where x = relative pressure (humidity, where water 
vapor is the adsorbate) expressed as a dimensionless 
decimal such that ^ x ^ 1. The dependent variable 
yis, in itself, an expression: 



y= x/[q (1 -x)] 



(2) 



containing both x and q. The latter quantity represents 
the number of moles of adsorbate per g of adsorbent 
obtained under steady state conditions at the corre- 
sponding value of x. 

The y-intercept, a, and slope, b 9 are the parameters 
ll(q m c) and (c — l)/(g m c), respectively. The quantity 
g m , in turn, represents the number of moles of ad- 
sorbate required to constitute a monolayer and, as 
such, is proportional to the surface area. The c-value 
is related exponentially to an average energy of 
adsorption. In the conventional application of eq (1), 
the parameters a and b are evaluated from the ex- 
perimentally obtained adsorption points (%, q) prefer- 
ably by use of the method of least squares [8]. These 
"best" estimates of a and b are then used to obtain 
q m and c. Finally, the surface area is obtained con- 
ventionally by multiplying q m by an appropriate 
constant sometimes referred to as the "packing 
factor." 

Most laboratory workers in surface chemistry are 
painfully aware of the fact that the good linearity of 
eq (1) is destroyed when the experimentally deter- 
mined initial weight of the sample is in error. One 
should, therefore, expect that the best straight line 
would result only when a zero-point weight very close 
to the correct value is used in the calculations. This 
criterion was employed in the present investigation 
for determining the dry weight of a solid adsorbent 
when water vapor was the adsorbate. In other words, 
the humidified sample weights themselves, utilizing 
a least-squares technique, were permitted to determine 
that value of the dry weight, p, along with the usual 
two parameters a and b required to optimize the 
linearity: 

If W i represents the weight (in grams) of the sample 
at the ith point corresponding to a relative humidity 
Xi, the number of moles of water vapor adsorbed per 
g of dry adsorbent would be 



18p ' 



(3) 



The value of x\ is assumed to be error-free here as it 
is in the conventional derivation. 

The residual, /?,-, is the experimental error of the 
ordinate of the £th point as indicated by the difference 



ji — {a + bxi). Equations 
combined to yield: 



(1), (2), and (3) may be 



Ri = 



I8x ip 



0-x i )(r»-p) 



— (a + bxi). 



(4) 



It is required to determine the values of the three 
parameters p, a, and b for which the sum of the 
squares of the residuals is a minimum. 2 

This amounts to solving three nonlinear simultaneous 
equations for these parameters given by eqs (5), (6), 
and (7): 



d Y n 
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—Ayr* =° 



(5) 
(6) 

(7) 



By taking the partial derivative of the summation 
with respect to p as indicated by eq (5), the result is 
eq (8): 

A *Kl lA y x ' Wi 

a ^ (i-*,k^-p) 2 £ (i-x t )(w,- P y 



n A Wi 

-^2 {1 - Xi y (Wi - p) s =0. (8) 



In the same way, eqs (6) and (7) give rise to eqs (9) and 
(10), respectively. 
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(9) 



(10) 



A cursory inspection of eqs (8), (9), and (10) discloses 
that the nonlinearity is ascribable entirely to the 
unknown parameter, p. In fact, the terms to the left 
of the equals signs in eqs (9) and (10), respectively, 
are identified with those obtained [8] in the "normal- 
ized" equations during the derivation of the least- 
square solution for a linear system. This observation 
suggests a relatively simple method for the simul- 
taneous solution of eqs (8), and (9), and (10) for the 
parameters p, a, and b. An initial value p\ is selected 
for p which is not quite as large as the lowest value of 
Wt. This value of p is substituted in eqs (9) and (10) 



2 This is rigorously valid only when the error represented by /?, has a normal (Gaussian) 
distribution over its entire range. Experience has confirmed that in the present application, 
this assumption is reasonably well approximated. 
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making it possible to solve them simultaneously for 
a and b which are then designated a\ and b\ % respec- 
tively. The values of a u bi, and p\ are then substituted 
in the left side of eq (8). The resultant value is called 
$1. 

Since p\ was initially selected to be somewhat high, 
it is decreased by a small predetermined amount to 
obtain its second iterative value, p 2 . Corresponding 
values of a 2 and ft 2 are again calculated from eqs (9) 
and (10), as before, and all three contribute to the 
determination of <$> 2 from eq (8). If <J> 2 has the same sign 
as 4>i, it should be smaller in absolute value and the 
iterative process is repeated until there is a change in 
sign of <£. At that point the incremental change in p is 
reduced to one-tenth of its previous value and its sign 
is automatically reversed. Again, the new direction is 
maintained with the smaller increment until the next 
change in sign of <f> is encountered. The process is 
repeated until negligible changes are obtained in suc- 
cessive values of p, a, and b while successive values 
of <£> approach zero in the limit. Usually this occurs 
within about 25 iterations and gives an estimate of 
surface area and goodness of fit as well as the dry 
weight of the sample. 

The entire process has been programmed for 
computer operation (FORTRAN IV) and is given in the 
appendix (fig. Al) in a form used with a teletype time- 
sharing terminal.' 1 

The values for the slope and the intercept were 
combined in the usual way at the end of each iteration 
by taking the reciprocal of their sum to obtain q,„. 
The factor 7.529 X 10 4 m 2 /mol was used [10] to convert 
q m to surface area. Only the final value, however, after 
convergence was complete could be considered valid. 

(2) The Harkins and Jura (H-J) isotherm equation. 

A linear relationship based on a two-dimensional 
equation of state was proposed by Harkins and 
Jura [3]. It took the form 



yhere 



B 



a> = ln (l/x) 



(id 

(12) 



while A and B are constants. 



There are no relative humidity restrictions on its 
range of applicability. (This is in contrast to the more 
familiar BET equation valid only within the inclusive 
values from about 0.1 to 0.3.) The surface area applica- 
ble to the H-J equation is proportional to the square 
root of the parameter A. The proportionality constant 
is, of course, dependent upon the adsorbate and the 



3 The Control Data "KRONOS" System (6000 series) high-speed digital computer re- 
quired only about 1 Vi s in a typical computation. The central computer time can be further 
shortened to about '/a s when several independent determinations are performed, since 
the compiling time (c.a. 1 s) need not be repeated. Similar results were obtained using the 
Sperry Rand "UNI VAC 1108" System. 

Certain commercial materials and equipment are identified in this paper in order to 
specify adequately the experimental procedure. In no case does such identification imply 
recommendation or endorsement by the National Bureau of Standards, nor does it imply 
that the material or equipment identified is necessarily the best available for this purpose. 



temperature of the isotherm. The value used in this 
work corresponding to water vapor at 23 °C was 
5.620 Xl0 4 m 2 /mol [12]. 

Using the same procedure described for the BET 
theory, the least square normalized equations re- 
sulting from the H-J isotherm equation are: 
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(13) 



(14) 



(W t -p) 
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Wi 



(Wt-p) 



(15) 



The iterative procedure previously described was 
again resorted to and a computer program given in 
appendix (fig. A2) was tailored to accommodate the 
H-J equation. This resulted in converged estimates 
of dry weight and surface area while at the same time 
affording a measure of goodness of fit. 

(3) Other Isotherm Equations 

The Polarization Equation of deBoer and Zwikkcr 
[9]: 



q = a +b' In w 
and the Polanyi Potential Theory [9]: 
q = a" + b"a>-V 3 



(16) 



(17) 



were both subjected to the same treatment as described 
for the other two theories and will be discussed in a 
later section. 

3. Equipment and Procedures 

The gravimetric adsorption equipment employed 
used a circulatory flow system and a means for 
supplying water vapor at predetermined humidities 
to the carrier gas as described in a recent publication 
[10]. A drying tube mounted vertically and filled with 
P 2 5 supported on glass wool had a trap sealed to its 
lower end for holding back the phosphoric acid pro- 
duced and isolating it from the gas stream during 
desiccation of a sample. The drying step was carried 
out at room temperature (23 °C) as were all of the 
humidity exposures. Drying was performed either ini- 
tially or at the end of a humidity sequence, depending 
upon whether it was an adsorption or desorption 
series. The zero-humidity data were, of course, not 
used in the calculations and served only as experi- 
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mental verification of the calculated values. Saturated 
salt solutions were used for obtaining the desired hu- 
midities, details of which were described previously 
[10]. However, any standard commercial equipment 
would work equally well in order to apply the tech- 
nique in any particular case. 

The order in which the experimental points are 
obtained is extremely important. They must be in a 
continually decreasing (desorption or drying) or in a 
continually increasing (adsorption) sequence to avoid 
hysteresis scanning loops [11]. 

Table 5 shows that even low area adsorbents 
respond well in dry weight prediction regardless of 
whether adsorption or desorption data are used. 

The dry weights predicted in table 6 are significant 
because most of the data were collected at relative 
humidities above 50 percent. The ability to do this is 
especially valuable in nondestructive testing where 
complete dehydration might cause irreversible changes 
in the sample. 

4. Results 

Tables 1 through 6 summarize the results of several 
determinations. The grouping of the experiments helps 
to emphasize the capabilities as well as the limita- 
tions inherent in the method. 

Most of the examples in table 1 give reasonably good 
agreement in surface area between BET and H-J. The 
experimentally measured dry weights, in the majority 



of cases, are slightly greater than the calculated values. 
This suggests the possibility that the last traces of 
physically adsorbed water may not have been com- 
pletely removed by the experimental procedure used. 
In one instance this effect is rather striking. 4 

The two experiments involving barium glass 5 
(water extracted) represent the initial and final adsorp- 
tion sequence in an adsorption-desorption-adsorption 
cycle in which the relative humidity started at zero — 
went stepwise to saturation — then down to zero— and 
finally back up again to 0.334. 

Good agreement between surface area predicted 
from desorption data with that from adsorption is 
not generally to be expected. Table 2 shows exceptions 
occur where hysteresis is either absent entirely or 
negligibly small. 

Some insight may be deduced from the results of 
both isotherm equations even when only desorption 
(drying) data are available. Fair agreement between 
the "areas" as shown in the latter two examples of 
table 3 where hysteresis is known to be virtually 
nonexistent suggests that these values may be reason- 
ably valid estimates of the true surface area. 

Table 4 is self explanatory and shows consistency 
of both surface areas and predicted dry weights 
between widely different portions of the adsorption 
isotherm. 



4 See Cellulose Acetate Membrane (table 4). 

5 Corning Glass Company X95-1FD. 



Table 1. Adsorption 

Low humidity range (0.082 -* 0.334) 
Comparison of BET with H-J for the Same Da 



Expt. No. 




Surface 
area 


Dry weight 




(m 2 /g) 


(g) 


A503094 


Synthetic Hydroxyapatite (6 points) 
Calculated: 

BET 


47.2 
44.4 


2708 




H-J 


2717 




Experimental 


.2733 


A444168 


Activated Bentonite (5 points) 
Calculated: 

BET 


181.4 
182.5 


.7161 




H-J 


.7193 




Experimental 


.7378 


A444161 


Polycarboxylate Cement (4 points) 
Calculated: 

BET 


60.4 
61.3 


1.255 




H-J 


1.241 




Experimental 


1.249 


A444143 


Zinc Phosphate Cement (6 points) 
Calculated: 

BET 


41.8 
32.4 


1.5646 




H-J 


1 5703 




Experimental 


1 5779 
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TABLE 1. Adsorption — Continued 



Expt. No. 




Surface 
area 


Dry weight 




(m 2 /g) 


(g) 


A444155 


Barium Glass — as received (6 points) 
Calculated: 

BET 


3.22 
3.15 


1.2886 




H-J 


1.2889 




Experimental 


1.2895 


A444156 


Barium Glass — Water Extracted (4 points) 
Calculated: 

BET 


8.7 
4.6 


0.9636 




H-J 


.9648 




Experimental 


.9649 


A444175 


Barium (Mass — Water Extracted (6 points) 
Calculated: 

BET (no convergence) 








H-J 


9.5 


.9646 






.9657 


A444164 


Zinc Oxide-Eugenol Cement (4 points) 








Calculated: 

BET (no convergence) 






H-J 


5.2 


3.1274 






3.1297 











FABLE 2. Adsorption versus desorpti 



Expt. No. 


Hysteresis: Absent or minimal 


Surface 
area 


Dry weight 




(m 2 /g) 


(g) 




Poly(methyl methacrylate) (Low Humidity range 0.082 —» 0.334) 






A444065 


Adsorption (5 points) 
Calculated: 








BET 


30.4 
24.9 


3.5029 




H-J 


3.4878 


A444064 


Desorption (6 points) 






Calculated: 








BET 


(31.6) 

(29.4) 


3.5022 




H-J 


3.4847 




Experimental 


3.5039 




"BIO GEL" (commercial hydroxyapatite) 






A444136 


Adsorption (Intermediate Humidity range 0.312-^ 0.804) 
Calculated (4 points): 








BET (not applicable) 


43.3 







H-J 


0.3788 


A444138 


Desorption (High Humidity range 0.926^0.542) 
Calculated (5 points): 








H-J 


(47.2) 


.3792 




Experimental 


.3802 
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Table 3. Desorption 

Low Humidity Range (0.334-^0.082) Comparison of BET with H-J for the same data 



Expt. No. 




Surface 
area 


Dry weight 




(m 2 /g) 


(g) 


A444159 


Barium Glass — Water Extracted (5 points) 
Calculated: 

BET 


(9.2) 
, (2.3) 


0.9639 




H-J 


.9659 






.9657 


A503137 


NPG-GMA Coated HA (6 points) 
Calculated: 

BET 


(14.6) 
(12.6) 


.4808 




H-J 


.4799 




Experimental 


.4805 


A444139 


"BIO GEL" (commercial hydroxyapatite) (5 points) 


(40.5) 
(59.4) 






Calculated: 

BET 


.3792 




H-J 


.3784 




Experimental 


.3811 











Table 4. Adsorption 

Low humidity versus high humidity 



Expt. No. 




Surface 
area 

(m 2 /g) 


Dry weight 

(g) 




Dentin 






A503140 


Calculated (6 points 0.082^0.334) 








BET 


146 


2.0202 




H-J (poor convergence) 








Experimental 


2.0671 


A503078 


Calculated (5 points 0.082^0.312) 








BET (no convergence) 


143 






H-J 


.2847 


A503079 


Calculated (4 points 0.334 -> 1.000) 






BET (not applicable) 


150 






H-J 


.2848 




Experimental 


.2898 




Cellulose Acetate Membrane 






A444148 


Calculated (6 points 0.082^0.334) 








BET 


191 
176 


.06424 




H-J 


.06480 


A444149 


Calculated (7 points 0.435^ 1.000) 










197 






H-J 


.06491 




Experimental 


.06681 











570 



Table 5. Adsorption — desorption 

High humidity range (0.435-0.971) 



Expt. No. 




Surface 
area 
(m 2 /g) 


Dry weight 
(g) 




Pulverized Dental Amalgam 






A444124 


Adsorption (7 points) 
Caleulated 








H-J 


0.44 


5.6135 


A444126 


Desorption (4 points) 






Calculated 








H-| 


(0.74) 


5.6132 






5.6145 

















TABLE 6. High humidify drying l desorption ) curves for predicting water content i 



Expt. No 



Surface 
area 
(m 2 /g) 



Dry weight 
(g) 



A444145 



A444153 



A444158 



A444171 



A503091 



A503136 



Cellulose Acetate Membrane (Humidity Range 0.926— > 0.655) 

Calculated (4 points): H-J 

Experimental 



Barium Class (Humidity Range 0.926 -> 0.435) 

Calculated (6 points): H-J 

Experimental 



Extracted Barium Class (Humidity Range 0.971 -* 0.435) 

Calculated (7 points): H-J 

Experimental 



Activated Bentonite (Humidity Range 0.926-* 0.542) 

Calculated (5 points): H-J 

Experimental 



Commercial Synthetic Hydroxyapatite (Humidity Range 1.000— > 0.757) 

Calculated (5 points): H-J 

Experimental 



NP(;-CMA Coated HA (Humidity Range 1.0O0-* 0.542) 

Calculated (7 points): H-J 

Experimental 



(148.3) 



(2.5) 



(3.4) 



(242) 



(117) 



(25) 



.06645 
.06681 



1.28897 
1.28947 



.9663 
.9657 



.7113 
.7378 



.2612 
.2733 



.4788 
.4805 



*Other examples included in tables 2 and 5. 



5. Discussion 

The fact that the least-squares technique can be 
applied successfully for determining unbound water 
content along with surface area for the BET and the 
H-J isotherm equations raises the question of whether 
it would work equally well for other adsorption theories. 
Accordingly, computer programs were developed for 
testing the Potential theory and the Polarization theory, 
each of which could be expressed in linear form [9]. 
It soon became clear that this method would not be 
applicable to either one of these latter isotherm 
equations. They seemed to be "ill-conditioned" or 
not sufficiently sensitive to perturbation in sample 
weight to be of practical use in instances where both 
BET and H-J converge readily. Further testing was 



therefore abandoned. The possibility remains that 
still other isotherm equations [2, 4, 5] (besides the four 
tested here) might well be adaptable to this technique. 
It is instructive to assess the relative advantages of 
BET and of H-J as applied to the present work: 

(1) The H-J would, of course, be used for all values 
beyond the relative pressure (about 0.3) of the appli- 
cable BET range. 

(2) For very small surface area samples where the 
points in the BET range may be too erratic for success- 
ful convergence, the points higher up on the isotherm 
may have sufficiently reduced relative error to give 
successful convergence using H-J. 

(3) In the valid BET region where both isotherm 
equations could be used, the BET equation appeared 
to be the more sensitive. When the experimental 
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error of the points was excessive, neither computer 
program would converge. At moderate error such that 
only one equation led to successful convergence it 
was usually H-J (see tables 1 and 4). 

(4) In many instances where no convergence was 
found when using the BET equation, a minimum in the 
residual variance (mean-square deviation of the or- 
dinates of the experimental points from the assumed 
curve) was useful as a rough estimate of the correct 
dry weight and surface area. 

(5) When convergence was attained using both 
isotherm equations and the predicted surface areas 
were in reasonably good agreement (see first four 
experiments in table 1), it served to reinforce one's 
confidence in the validity of the results. When agree- 
ment was somewhat less than desired for a homo- 
geneous adsorbent (see table 1, Expt A444156), the 
BET value would be preferred. Note that the BET 
surface area value of 8.7 m 2 /g (rather than the H-J 
value of 4.6) was in better agreement with the results 
of Experiment A444175 for the same adsorbent (also 
table 1) for which the H-J value of 9.5 m 2 /g was 
calculated. 

(6) For a nonhomogeneous adsorbent where aggregate 
behavior might lead to an erroneous BET value [12] as 
a result of large differences in energies of adsorption 
among the components of the solid, the H-J area would 
always be preferable. 

Throughout tables 2, 3, 5, and 6 surface area 
values obtained from desorption experiments are shown 
in parenthesis. They are likely to be erroneous unless 
the desorption curve is substantially a retracing of the 
adsorption curve (absence of hysteresis). This is a 
property of the adsorbent-adsorbate system and must 
be considered case by case. Virtual absence of 
hysteresis is portrayed by the materials listed in table 
2 and exists to only a slight degree for those listed in 
table 3. The dry weight values, however, are just as 
valid with desorption as with adsorption provided the 
desorption originates at saturation resulting in a 
descending boundary curve [13]. Other instances which 
have been found to agree with the experimental re- 
sults for dry weight include primary descending 
scanning curves; although the predicted surface area 
would not be valid. A case in point is illustrated by an 
experiment not shown in the tables. Water vapor was 
adsorbed on a desiccated sample of Bio Gel 6 to a 
relative humidity of 0.804 after which seven additional 
successive desorption values were obtained. Con- 
vergence was successful using H-J resulting in a 
predicted dry weight of 0.3788 g compared with the 
experimental value of 0.3811 g obtained after taking 
the sample again to dryness. 

Desorption along scanning curves other than those 
already mentioned as well as adsorption along scanning 
curves, will invariably result in erroneous values for 
dry weight prediction in spite of satisfactory con- 
vergence of the computer program. Two experiments 
performed with the same sample of human dentin will 
serve to illustrate this. In both cases, the starting 



B BIO-RAD Laboratories, Richmond, California. 



point was saturation (relative humidity = 1.00). In the 
first experiment (A444090), desorption occurred along 
the descending boundary curve to a value of 0.175 
relative humidity. This was followed by adsorption 
along a primary ascending scanning curve to a rela- 
tive humidity of 0.757. Desorption was next carried 
out stepwise along this secondary descending scan- 
ning curve. The sample weights measured were: 
0.249303 g, 0.247758 g, 0.246568 g, and 0.245402 g at 
the corresponding relative humidities of 0.757, 0.655, 
0.542, and 0.435, respectively. The dry weight pre- 
dicted from using the H-J equation was 0.2356 g as 
compared with 0.2288 g measured experimentally. In 
the other experiment (A444091) the initial turning 
point between the descending boundary curve and 
the primary ascending scanning curve took place at 
0.122 r.h. (instead of 0.175 r.h.). In this experiment, 
the sample weights at the same respective humidities 
along the secondary descending scanning curve were 
0.248456 g, 0.247364 g, 0.245881 g, and 0.244989 g. 
These values gave rise to a predicted dry weight of 
0.2369 g. It is significant that in both of these experi- 
ments the dry weight was predicted about four percent 
on the high side compared with the actual measure- 
ment. A detailed study of the behavior of scanning 
curves in adsorption-desorption hysteresis is beyond 
the scope of this paper [13]. 

It was stated earlier in this discussion (subparagraph 
3) that the BET equation was extremely sensitive 
regarding the effect of experimental error on con- 
vergence of the computer program. If a sufficient 
number of points are available in a system in which 
the BET equation is known to be valid and yet the 
computer program does not converge, it is useful to 
plot the original data especially when it is impractical 
(or impossible) to repeat the experiment. If one of the 
points does not fit a smooth curve which passes through 
the remaining points, those remaining points may very 
well allow the computer program to converge. Of 
course, the results should be examined carefully for 
reasonableness of the predicted quantities before they 
are accepted. This suggestion is submitted in lieu of 
a valid test for outliers [14] in a system as compli- 
cated as this. 

It is appropriate to caution the reader that the 
method described in this paper is not proposed as a 
substitute for measuring the zero humidity point in 
routine surface area determinations. It is also en- 
tirely feasible that experimental errors might combine 
in such a way as to permit convergence and yet yield 
incorrect results. In other words, the need for careful 
experimental work is probably even more important 
when the additional demand is placed on the sorption 
points to estimate the zero humidity value as well as 
the surface area. 

Some applications in which the methods described 
in this paper would be especially valuable are: 

(1) Where routine measurements have yielded good 
linear BET results yet, now and then, a sample of the 
same type of material when plotted in accordance 
with eq (1) produced a curve instead of a straight line. 
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In such instances, the entire experiment would 
normally have to be discarded and all because an 
error in the experimental determination of the dry 
weight was reflected in all the other points on the plot. 

(2) Dry weight of materials that evolve gases during 
decomposition at low pressures. 

(3) In agricultural products (such as grain) and 
other materials where the last traces of water may be 
extremely slow coming off during drying (desorption), 
provided again that the isotherm equation employed 
is known to be appropriate for that system. 

(4) In biological applications where viability of the 
organism (or tissue) would be destroyed by complete 
drying while, perhaps, equilibration at, say, fifty 
percent relative humidity might be tolerated, non- 
destructive determinations of dry weight could be 
made. This would be especially useful in growth rate 
studies. 

(5) In situations where no other method may be 
available for estimation of dry weight and/or surface 
area from previously obtained data which may be 
incomplete. 

6. Summary 

(1) It has been demonstrated that adsorption or 
desorption data are capable of estimating the unbound 
water content of solid adsorbents. This was accomp- 
lished by the convergence of a computer program 
which employed a least-squares technique applied 
to an appropriate isotherm equation. 

(2) The Brunauer, Emmett, and Teller free surface 
equation and the Harkins and Jura equation were 
readily adaptable to this treatment (the Polarization 
equation and the Potential theory equation did not 
respond). 

(3) A valid estimate of surface area was also obtained 
when data from adsorption experiments were used or 
with such desorption experiments for which no 
hysteresis existed. Where there was hysteresis, the 
predicted surface area while incorrect was a useful 
order-of-magnitude estimation, usually as an upper 
limit. 

(4) Applications of the technique for determining 
the dry weights were compared with experimentally 
measured values for solids which cover a wide range 
in surface area. The agreement was within about one 
percent and the calculated dry weight was usually on 
the low side. 

(5) The capabilities as well as the limitation in the 
application of these methods were compared and 
discussed. 

7. Appendix. Procedure for Using Isotherm 
Data With Time Sharing Teletype Terminal 



Jura equation. A minimum of four data points is 
required, although a larger number is preferred. The 
adsorption (or desorption) data can be pre-punched 
on paper tape before connecting the terminal with the 
computer. (The input format has been designed such 
that the same data tape may be used with either 
program.) A typical data file, for example, (A444153) 
transposed from the tape is shown in table Al. The 
initial line (here designated as 10) contains: 

a. The number of data points in the experiment 

(in this case, 6). 

b. The maximum number of iterations. 

c. A starting value (tentative) for the dry weight. 

d. The initial decrement (or amount by which the 

preceding numeral is changed in initiating the 
iterative process). 

Subsequent lines of table Al contain the sample 
weight along with the corresponding humidity to 
which the sample was exposed until that steady-state 
weight had been attained. 

The order of the points in table Al indicates in this 
illustration that a desorption process originating at a 
relative humidity of 1.00 was measured initially at 
0.926 r.h., and finally at 0.435. While a maximum of 
33 iterations was provided for in this example, only 
about 20 were needed before no further change of any 
consequence resulted. 

Table Al. Input data for Expt. A444153 

10 06 33 01.289000-0.000400 

11 01.294267 0.926 

12 01.291396 0.804 

13 01.290930 0.757 

14 01.290616 0.655 

15 01.290353 0.542 

16 01.290112 0.435 

The initiating value of 01.289000 g was chosen to be 
less than any of the data points and yet obviously 
greater than the expected dry weight. The choice of 
—0.000400 for the weight change meant that a value of 
1.288600 g for a tentative dry weight would be tested 
by the computer in its next iteration. The Harkins and 
Jura program (HARJUR) would, of course, be used for 
the data file A444153, since the humidities are beyond 
the valid range applicable to the BET theory. 

In the course of the computer run the value of 4> 
must change signs repeatedly as it converges toward 
zero. This is a necessary criterion for a successful 
determination. The value predicted for dry weight 
for this particular set of data was 1.28897 g while the 
predicted 7 surface area was 2.49 m 2 /g. (Other columns 
present in the computer printout (not shown here) 
represent intermediate steps in the calculations or 
estimate the goodness of fit to the linear form of the 
isotherm equation.) 

An interesting application of the computer program 
is the degenerate case where the iterations are sup- 
pressed. This would occur when the dry weight is 



Kortran ComOUter DrO°TaniS are ffiven in fiffUre Al 7 The predicted surface area would normally not be valid in the case of a desorption ex- 

,, . . .5? A r» r l tt l • l periment (such as this) unless hysteresis was absent. It is useful, however, as an estimate of 
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00010 program vasser c input, output* taped 

00020*this program i£ designed to find the dry weight pcjj) vheke 
00030*xci) is the relative humidity, vci) is the veight of sample 
00040*at equilierum vi th xci), and j is the iteration of the par- 

000 50*ameter. max i is ii. max j is jj. the value of ii 
00060*mvst ee at least 4. the range of x(i) should ee 0.08 to 0.33. 
00070 dimension vc50), xc50), pc99), f3c50), f4c50), f5c50), 

00080+ F6C50), F8(50) 

00090 READ CI, 61) II, JJ, PCI), DELP 

00100 61 FORMAT C3X, 12, IX, 12, IX, F9.6, IX, F9.6) 

001 10 PFINT 62 

00120 62 FORMAT C/3X, 2HII, IX, 2HJJ, 5X, 4HPC1>, 7X, 4HDELP//) 

00130 PRINT 61, II, JJ, PCI), DELP 

00140 READ CI, 63) CVCI), XCI), 1=1,11) 

00150 63 FORMAT C 3X, F9.6, IX, F5. 3) 

00160 PRINT 64 

00170 64 FORMAT C///5X, 4HVCI), 5X, 4HXCI)//) 

00180 PRINT 63, CVCI), XCI), 1=1,11) 

00190 PRINT 65 

00200 65 FORMAT C////1X, 1HJ, 6X, 4HPCJ), 8X, 3HPHI, 9X, 1HA, 9X, 

00210+ 1HE, 8X, 4HAREA, 7X, 1 0HMN SGR DEV//) 

00220 All = II 

00230 PPEPHI =0.0 

00240 DO 104 J=1,JJ 

00250 GX2 =0.0 

00260 DO 92 1=1,11 

00270 92 GX2 = CX2 + XCI)**2 

00280 GX = 0.0 

00290 DO 94 1= 1,11 

00300 94 GX = GX + XCI ) 

00310 DEN = AII+GX2 - GX**2 

00320 DO 111 1=1, II 

00330 111 F3CI) = XCD/C 1 . 0-X C I ) ) / C V C I ) -PC J) ) 

00340 GF3 = 0. 

00350 DO 112 1=1,11 

00360 112 GF3 = GF3 + F3CI) 

00370 DO 1 13 1=1, II 

00380 113 F4CI) = XCI)*F3CI) 

00390 GF4 = 0.0 

00400 DO 114 1=1,11 

00410 114 GF4 = GF4 ♦ F4CI) 

00420 DO 115 1=1,11 

00430 115 F5CI) = F3CI)*VCI)/CVCI)-PCJ)) 

00440 GF5 =0.0 

00450 DO 116 1=1, II 

00460 116 GF5 = GF5 + F5CI) 

00470 DO 117 1=1,11 

00480 117 F6CI) = XCI)*F5CI) 

00490 GF6 =0.0 

00500 DO 118 1=1,11 

00510 118GF6 = GF6 + F6CI) 

00520 DO 1 19 1=1,11 

00530 119 F8CI) = F3CI)*F5CI) 

00540 GF8 = 0.0 

00550 DO 120 1=1,11 

00560 120 GF8 = GF8 + F8 C I ) 

00570 A = 18.0*PCJ)*CGX2*GF3 - GX*GF4)/DEN 

00580 E = 18.0*PCJ)*CAII*GF4 - GX*GF3)/DEN 

00590 AREA = 7. 529E+04/ C A+E) 

00600 PHI = A*GF5 + B*GF6 - 1 8 . 0*PC J) *GF8 

00610 SSQRS = 0.0 

00620 DO 4 1=1,11 

00630 4 SSGRS = SSQRS + CVCI) - PCJ)*C1.0 + 18.0*XCI)/ 

00640+ C1.0 - XCI))/CA + E*XCI))))**2 

00650* SSQRS OF DEVS OF VCI) FROM CURVELINEAR REGRESSION 

00660 VAR = SSCRS/CAII - 3.0) 

00670 PRINT 99, J, PC J), PHI, A, E, AREA, VAR 

00680 99 FORMAT C 1 2, 2EI3.5, 2E10.2, 2E12.4) 

00690 IF CPHI*PEEPHI) 95,96,96 

00700 95 DELP = -0.1*DELP 

00710 96 PREPHI = PHI 

00720 PCJ+1) = PCJ) + DELP 

00730 104 CONTINUE 

00740 STOP 

00750 END 

FIGURE Al. Computer Program, WASSER, (FORTRAN IV) for determination of unbound water content (dry weight) 
and Surface Area from water vapor sorption data according to the BET adsorption theory. 
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000 10 PROGRAM HAFJUF.C INPUT, OUTPUT, TAPED 

00020+THIS PROGRAM IS DESIGNED TO FIND THE DRY VEIGHT P(JJ> 
00030*VHERE X(I) IS THE RELATIVE HUMIDITY, V(I) IS THE 
00e40*' T EIGHT OF SAMPLE AT ECUILIERUM VITH XCI), AND J IS THE 
00050*ITERATION OF THE PARAMETER* MAX I IS II. MAX J IS JJ. 
00060*THE VALUE OF II MUST EE AT LEAST 4. 

00070 DIMENSION VC50), X(50), PC99), H1C50), H2C50) 

00080 READ(1,6D II, JJ, PCD, DELP 

00090 61 FORMAT (3X, I 2, IX, 12, IX, F9.6, IX, F9.6) 

00100 PRINT 62 

00110 62 FORMAT ( /3X, 2HI I , IX, 2HJJ, 5X, 4HFC 1 ) , 7X, 4HDELP// ) 

00120 PRINT 61,- II, JJ, PCD, DELP 

00130 EEADCD63) CVCD,XCI), 1=1,11) 

00140 63 FOFMAT C 3X, R9 . 6, IX, R5. 3 ) 

00150 PRINT 64 

00160 64 FORMAT C///5X, 4HVCI), 5X, 4HXCI)//) 

00170 PRINT 63, CVCI), XCI), 1=1,11) 

00180 PRINT 65 

00190 65 ROFMAT C////1X, 1HJ, 6X, 4HPCJ), 8X, 3HPHI, 9X, 
00200+ 1HA, 9X, 1HE, 8X, 4HAFEA, 7X, 1 0HMN SCR DEV///) 

00210 All = II 

00220 PREPHI =0.0 

00230 DO 104 J = 1,JJ 

00240 DO 90 I = 1, II 

00250 90 H1CI) = ALOGC 1. 0/XCI ) ) 

00260 GH1 = H1C 1 ) 

00270 DO 91 I = 2, II 

00280 91 GH1 * GH1 ♦ H1CI) 

00290 DO 92 I = 1, II 

00300 92 H2CI) = 1.0/CVCI) - PCJ>) 

00310 GH2 = H2C 1 )**2 

00320 DO 93 I = 2,11 

00330 93 GH2 = GH2 + H2CI>**2 

00340 GH12 = HI C D*H2C 1 )**2 

00350 CO 94 I = 2,11 

00360 94 GH12 = GH 1 2 ♦ H 1 C I ) *H2C I ) **2 

00370 GH4 = H2C 1 )**4 

00380 DO 95 I = 2,11 

00390 95 GH4 = GH4 + H2CI)**4 

00400 GHV3 = VC 1 )*H2C 1 )**3 

00410 DO 96 I = 2,11 

00420 96 GHV3 = GHV3 ♦ VC I ) *H2C I ) **3 

00430 GH1V3 = VC D*H1C D*H2C 1)**3 

00440 DO 97 I = 2, 1 1 

00450 97 GH1V3 = GH1V3 + VC I ) *H 1 C I ) *H2C I ) **3 

00460 GHV5 = VC 1 )*H2C 1 )**5 

00470 DO 98 I = 2, 1 1 

00480 98 GHV5 = 0HV5 ♦ VC I ) *H2C I ) ** 5 

00490 DEN = C 18. 0*PC J) )**2*CAI I*GH4 - GH2**2> 

00500 ANUM = AII*GH12 - GH1*GH2 

00510 A = ANUM /DEN 

00520 ENUM = C 1 8 . 0*PC J ) ) **2*C GH2*GH 1 2 - GH1*GH4) 

00530 E = ENUM/DEN 

00540 PHI = A*PC J)*CGH1V3 - A* C 18. 0*PC J) ) **2*GHV5 •»■ E*GHV3) 

00550 SSQRS =0.0 

00560 DO 4 I = 1,11 

00570 DISC = A/CH1CI) + E) 

00580 IF CDISC) 70,69,69 

00590 69 DSCR = DI SC 

00600 SSQRS = SSQRS + C VC I ) -PC J)*C 1 . 0+ 18. 0*SGRTC DSCR) ) ) **2 

00610 GO TO 4 

00620 70 SSQRS = 1.0E+08 

00630 4 CONTINUE 

00640* SSQRS OF DEVS OF VCI) FROMM CURVILINEAR REGRESSION 

00650 VAP = SSQRS/CAII-3.0) 

00660 IF CA.LT.0.0) GO TO 111 

00670 AREA = 5. 62E+04* SQRTC A) 

00680 GO TO 112 

0069 1 1 1 AREA =0.0 

00700 112 PRINT 99, J, PCJ), PHI, A, E, AREA, VAR 

00710 99 FORMAT C I 2, 2E13.5, 2E10.2, 2E12.4) 

00720 IF CJ.EQ. 1) GO TO 101 

00730 IF CPHI*PREPHI) 100,100,101 

00740 100 DELP = -0. 1*DELP 

00750 101 PREPHI = PHI 

00760 PCJ+D = PCJ) + DELP 

007 7 104 CONTINUE 

00780 STOP 

00790 END 
Figure A2. Computer Program, HARJUR, (FORTRAN IV) for determination of unbound water content {dry weight) 
and Surface Area from water vapor sorption data according to the H-J adsorption theory. 
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available experimentally along with the rest of the 
data points and the experimenter is primarily interested 
in the least-squares value of the surface area. In the 
first line of the input data file (table Al) he would set 
the number of iterations equal to 01 and would insert 
the known dry weight for the initial value of p which 
is the next entry on that same first line. The final entry 
(corresponding to the initial decrement) while it will 
be ignored by the computer should, never-the-less, be 
present even if all its digits are zeroes. 

The computer programs used in this work are 
reproduced in figures Al and A2. 
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